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ABSTRACT 

Establishing the origin of the short-lived radionuclide (SLR) 26 Al, which was present in refractory 
inclusions in primitive meteorites, has profound implications for the astrophysical context of solar 
system formation. Recent observations that 26 Al was homogeneously distributed in the inner solar 
system prove that this SLR has a stellar origin. In this Letter, we address the issue of the incorporation 
of hot 26 Al-rich stellar ejecta into the cold protosolar nebula. We first show that the 26 Al atoms 
produced by a population of massive stars in an OB association cannot be injected into protostellar 
cores with enough efficiency. We then show that this SLR likely originated in a Wolf-Rayet star that 
escaped from its parent cluster and interacted with a neighboring molecular cloud. The explosion of 
this runaway star as a supernova probably triggered the formation of the solar system. This scenario 
also accounts for the meteoritic abundance of 41 Ca. 

Subject headings: ISM: bubbles — nuclear reactions, nucleosynthesis, abundances — stars: formation 
— stars: Wolf-Rayet 



1. INTRODUCTION 

Mor e than 30 years after the discovery of ILee et afl 
(|1976f ) that calcium-aluminum-rich inclusions (CAIs) 
from the Allende meteorite contained 26 Al (mean lifetime 
t-26 = 1-03 x 10 6 yr), the origin of this short-lived radionu- 
clide (SLR) remains an open question. High precision 
Mg isotopic analy ses of asteroids and bulk rocks from 
terrestrial planets (|Thrane et al.ll2006t ), as well as recent 
micrometer-scale data in chondrules (|Villeneuve et al.l 
2009), showed that 26 Al was homogenously distributed 
over at least the inner part of the solar system, i.e., over a 
reservoir of mass > 2 , and that no significant amount 
of freshly made 26 Al was added to the protoplanetary 
disk after the CAI formation. The maximum amount of 
26 Al that could have been synthesized by in situ parti- 
cle irradiation during the short dur ation of CAI forma- 
tion (~ 10 5 yr; IBizzarro et"aTll2004[ ) can account for the 
canonical 26 A1/ 27 A1 = 5 x 10~ 5 over a roc ky reservoir 
of only -0.1 M e (|Duprat fc Tatischeffl[2007l ). Thus, the 
origin of this SLR cannot be related to the nonthermal 
activity of the young Sun and has to be searched for in a 
stellar nucleosynthetic event contemporary with the for- 
mation of the solar system. 

An origin of SLRs in an a symptotic giant branch 
(AGB) star has been propo sed (Wasscrb urg et al.l 119941 ; 
ITrigo- Rodriguez et al.1 12009). but AGB stars are not as- 
sociated with star-forming regions and the probability 
of a chance encounter between an AG B star and a star- 
formi ng molecular cloud is very low ([Kastner fc Mversl 
[1991 . It is more likely that the protosolar system was 
contaminated by material freshly ejected from a mas 



sive s tar, either a Type II SN (e.g. Cameron fe Truran 
19771) or a Wolf-Rayet (WR) star (e.g. lArnould et al 



19971) . Massive stars have a profound influence on the 



surrounding m olecular clouds an d the process of star 
formation (e.g. ILee fe Chenll2007D . iCameron fe Truranl 



(|1977f) first suggested that a supernova (SN) responsi- 
ble for injecting SLRs into the presolar nebula may also 
have been responsible for triggering the formation of 
the solar system. Detailed numerical simulations have 
shown that such simultaneous triggering and injection 
is possi ble, but the inje ction efficiency is lower than re- 
quired (|Boss et al.ll2010fl . Alternatively, it has been sug- 
gested that a nearby SN (~0.3 pc) may have injected 
SLRs into the alre ady-formed protoplan etary disk of the 
solar system (see lOuellette et al.l I2007L and references 
therein). In this scenario, it is assumed that the Sun was 
born in a large stellar cluster containing massive stars. 
But this model is questionable, because (1) protoplan- 
etary disks in the vicinity of massive stars are exposed 
to a rapid photoevaporation and (2) the main-sequence 
lifetime of even the most massive stars is too long as 
compared to the mean life time of protoplanetary disks 
dGounelle fe Meibqmll2q08f) . 

iGounelle et al.l (|2009l i and iGaidos et all (|2009l ) re- 
cently suggested that the Sun is born in a stellar clus- 
ter of second generation, whose formation was triggered 
by the activity of a neighboring OB association. There 
are many observations of OB associations divide d in spa- 
tially separated subgroups of different ages (e.g. iBlaauwl 
1964), as well as observations of young stel lar objects 
locat ed on the border of H II regions (e.g. | Karr et all 
120091 ) and superbubbles (e.g. ILee fc Chenll2009f ) AdoptT 
ing such an astrophysical context, we study in this Letter 
how hot stellar debris enriched in 26 Al could be injected 
into a cold protostellar nebula. We show in Section 2 
that 26 Al produced by a population of massive stars in 
an OB association may not be delivered into molecular 
cores efficiently enough. We then study in Section 3 a 
possibility already m entioned in the pione ering work of 
lArnould et al.l (|1997t ) and more recently by IGaidos et al.l 
(2009) that the presolar nebula was contaminated by 
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Fig. 1. — Sketch illustrating the ejection of 28 Al by (a) a pop- 
ulation of massive stars in an OB association and (b) a WR star 
"running away" from its parent cluster. 

26 Al produced by a WR star that escaped from its parent 
cluster. 

2. 26 AL PRODUCTION BY AN OB ASSOCIATION IN A 
SUPERBUBBLE 

Most massive stars are born in OB associations, where 
multiple stellar winds merge and expand to form large 
cavities of hot ga s known as superbubbles (see, e.g., 
IParizot et al J 12004). The subsequent SNe generally ex- 
plode inside the wind-generated superbubble. The radius 
of a superbubble ca n be estimated fro m the standard 
wind bubble theory ((Weaver et al.lll977l ): 
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where tj&yi is the time in units of Myr after the on- 
set of massive star formation (assumed to be coeval 
for all stars), N* t aa = iV*/30 where N* is the num- 
ber of massive stars in the 8-120 Mq mass range, and 
n H,ioo = n n/ (100 cm' 3 ) where jih is the mean H number 
density in the ambient interstellar medium. The super- 
bubble radius is generally given as a function of the stel- 
lar wind mechanical power, L,„, instead of the n umber 
of massive stars (e.g. iMac Low fe; McCravl Il988[ ). But 
Equation (P) uses the recent result of IVoss et al.1 ()2009l ) 
that the mean wind power per star from a coeval pop- 
ulation of massive stars is nearly constant with time for 
~5 Myr and amounts to wl.5 x 10 36 erg s _1 . Simi- 
larly, the characteristic temperature and H number den- 
sity in the interior of a superbubble can be w ritten as 
(IWeaver et al.lfl977t [Mac Low fc McCravlll988tt 
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WR wind and SN ejections of Al occur at ^Myr ~ 3 
(jVoss et al.l 120091 ) . when the superbubble blown by the 
winds from the main-sequence stars has already reached 
a radius of several tens of pc (Equation (p}). Noteworthy, 
SN blast waves within a superbubble will usually become 
subsonic in the hot gas before they reach the supershell 
of swept- up interstellar material (|Mac Low fc McCravl 
1988; Parizot et aTll2004f ). This is true as well for winds 
of WR stars. Thus, most nuclei synthesized in mas- 
sive stars first thermalize in the hot superbubble inte- 
rior. Further incorporation of this material into molec- 
ular clo uds and star-forming s ystems takes more than 
10 Myr (jMever fc Clavtonll2000D . by which time the 26 Al 
will have decayed. 



To solve this issue, IGaidos et al.1 (|2009l ) proposed that 
26 Al ejected in WR winds can be rapidly incorporated 
into high-speed (^1000 km s -1 ) refractory dust grains 
of ^0.01-0.1 /zm size, that could dynamically decouple 
from the shocked wind gas and imbed themselves into 
the surrounding molecular material. But this proposal 
has two shortcomings. First, WR stars are thought to 
be a major contributor to the Galactic 26 Al detected 
through its gamma-ray decay line at E~ t — 1809 keV, and 
high-resolution spectroscopic observations of this emis- 
sion with the RHESSI and INTEGRAL gamma-ray satel- 
lites have shown that the line is narrow, A_B 7 = 1 
2 keV FW HM, consistent wi th the instrumental res- 
olution (see iDiehl et al.l [20061 and references therein). 
The non-detection of Doppler broadening in the Galactic 
1809-keV line provides an upper limit on the mean veloc- 
ity of the emitting 26 Al nuclei: w max ~ O.bcAE^/E^ ~ 
150 km s _1 (here, c is the speed of light). This max- 
imum velocity is much lower than the speed that dust 
grains must acquire to survive spu ttering as they pas s 
the WR wind termination shock (jGaidos et al.1 [2009) . 
Secondly, most grains formed in WR winds will slow 
down and stop in the superbubble interior before reach- 
in g the supershe ll. According to the classical estimate 
of ISpitzerl (laZg ), the range of a grain of size a gr and 
typical density p gr ~ 2 g cm~ 3 is X gr = a gr p gr = 
(2 x 10~ 6 g cm~ 2 )(a gr /0.01 /im). In comparison, the 
radial path length in a superbubble is 



X«B = 1.4?m 



n(r)dr 



(4.6 x 10~ 5 g cm" 2 
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where ?»h is the H mass and n(r) = Usb[1— (t/Rsb)] 
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(Weaver et al.l I1977T) . Thus, grains with 



0.2 jim 



do not reach the supershell. In fact, even much larger 
grains should stop in the superbubble interior, because 
the Spitzer formula can largely overesti mate the ran ge of 
interstellar dust grains in hot plasmas ([Ragotll2002T ). 

Dense clumps of molecular gas can be engulfed by the 
growing superbubble, if th ey were not swept u p by the 
expanding supershell (e.g. IParizot e t al. 2004). These 
clumps could potentially be enriched in 26 Al synthe- 
sized by WR stars and Type II SNe in the OB associa- 
tion. But recent two-dimensio nal hydrodynamic simula- 
tions (Bo ss et all I2008L 120101 ) suggest that the amount 
of 26 Al that could be injected into such a molecular 
cloud core would be too low to explain the solar sys- 
tem's canonical 26 A1/ 27 A1 ratio. Boss et al. found that 
only 2-5 xlO -5 Mq of hot SN shock front material could 
be incorporated into a cold molecular clump. But a 
1 Mq presolar cloud would need to be contaminated 
by ~10~ 4 M Q of SN matter to explai n the 26 Al mete- 
oritic abundance (Takigawa et al.ll2008f ). Although these 
two estimates are close, the main issue lies in the short 
lifetime of a small molecular cloud embedded in a hot 
plasma: the lifetime of a 1-M Q cloud against evapora- 
tion in the >10 6 K (Equation |2j)) superbubble interior is 
only ~10 5 yr (|McKee & Cowid[T977T) . much shorter than 
the duration of stellar main sequence. This scenario is 
therefore highly improbable. 



A runaway WR star as the origin of 26 Al in the early solar system 



3 



3. 26 AL PRODUCTION BY A RUNAWAY WR STAR 

If the vast majority, if not all O-type stars (the 
main-seq uence progenitors of WR stars) form in clus- 
ters (e.g. lLada fe Ladal 120031 ) . nearly half of them ac- 
quire velocities exceedin g the escape velocity from the 
cluster's potential well (jStond 119911 ). These runaway 
star£3 can be accelerated either by dynamical interac- 
tions wlJli_other_starsJii thedense cores of young clus- 
ters ((Leonard fc Duncan! I1990D or by the SN e xplosion 
of a companion star in a massive binary system (|Blaauwl 
fl96lh . A star moving with a velocity V* ^ 15 km s 
relative to its parent cluster leaves the associated super- 
bubble in less than 3 Myr (see Eq. [1]). About 20% 
of the O-ty pe stars have pecu liar velocities exceeding 
15 km s" 1 |de Wit et all 120051 ). These runaway short- 
living stars may have a significant probability of inter- 
acting with their parent molecular cloud complex. Out- 
side the hot gas, the star's motion is supersonic with 
respect to the ambient mediu m, which generates a bow 
shock (Ivan Buren et all [19901. There are many obser- 
vations of bow shocks created by runaway OB stars 
in the vicinity of young clus ters and associations (e.g. 
iGvaramadze fe Bom ans 20Ql). 

The form of a bow shock is determined by the balance 
between the ram pressure of the stellar wind and the 
ram pressure of the ongoing circumstellar (CS) gas. The 
pressure equilibrium is reached in the star's direction of 
motion at the so-called standoff distance from the star 
(|van Buren et al.lll~99Ct see Figure 1(b)), 



i?o = (0.92 pc)M, 
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where Mw—5 is the stellar wind mass-loss rate in units 
of 10~ 5 Mq yr" 1 , V w = W.1500 X 1500 km s" 1 is the 
wind's terminal velocity, and V*, 2 o = V*/(20 km s _1 ). 
The interaction region between the stellar wind and the 
CS medium is a shell bounded by two shocks, in which 
the flows slow down from supersonic to subsonic veloci- 
ties. One verifies that for Vw 3> V* the shell's mass is 
mainly due to the shocked CS gas. The contact surface 
between the shocked stellar wind and CS gas is unstable 
due to both the Kelvin-Helmholtz and R ayleigh- Taylor 
instabilities (Brig henti fe D'Ercoldll9 95b). such that we 
expect an efficient mixing of the wind-ejected material 
with the swept-up CS gas. 

The 26 Al/ 2 ' Al ratio in the bow shock shell of a runaway 
WR star just prior to the SN explosion can be estimated 
as a function of the polar angle 9 from the star's direction 
of motion (see Figure 1(b)) by 
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where A 2 6 is the total number of 26 Al nuclei ejected in 
the WR wind, / 26 = (t 26 /A W r) X (1 - exp(-A W R/T 26 )) 
is a factor that takes into account the decay of 26 Al dur- 
ing the duration Awr of the WR phase (Awr — 0.2- 

1 While traditionally the minimum pec uliar v elocity for classi- 
fying a star as a "runaway" is 40 km s — 1 {Blaauw 1961]), here we 
use this ter m for any st ar that has escaped from its parent clus- 
ter (see alsolStonc 1991]). The escape velocity ranges from several 
km s — 1 for loose and low-mass clusters to several tens of km s — 1 
for compact and massive ones. 
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Fig. 2.— Mean 26 A1/ 27 A1 ratio (thick curves, left axis) and H 
column density (thin curves, right axis) in the bow shock shell of 
a runaway WR star, as a function of position 9 in the shell (see 
Figure 1(b)), for V* = 20 km s" 1 and n H = 100 cm" 3 . The 
26 A1/ 27 A1 ratio is not shown for the 40 Mq star (see the text); 
besides, the Eg curve for this star merges with that for the 60 and 
85 Mq stars. 

1.4 Myr depending on the star's mass), X27 = 3.46 x 10~ 6 
is the 27 A1 abundance by number in the CS med ium as- 
sumed to be of solar composition (|Loddersll2003D . Rs{&) 
is the shell's radius (Figure 1(b)), and Tis(Q) is the 
shell' s H column density (£<?(0°) ~ 0.75-Ron H ; see lWilkinl 
Il996f) . W e took the 26 Al yields from the rotating stellar 
models of iPalacios etH] p005h . R s (0) and Z s (6) were 
calculated from the analytic solutions found by Wilkin 
(1996) in the thin-shell approximation. The WR star 

parameters Mw, Vw, and Awr were extracted from 
the grids of ro tating stellar mod els of Meynet & Maeder 
(2003; see also lVoss et"aTll2009D . 

Calculated 26 A1/ 27 A1 ratios in bow shock shells of run- 
away WR stars are shown in Figure 2 for stars of ini- 
tial masses 25, 60, and 85 M B (the 26 Al yield for the 
40 M Q star in not listed in IPalacios et alj|2005l) . The 
26 A1/ 27 A1 ratio is only weakly dependent on the star's 
initial mass, because both A 2 6 and the amount of swept- 
up CS matter (oc i?|Es) increase with increasing stellar 
mass. We see that ( 26 A1/ 27 A1) S reaches -1-2 xlO~ 2 at 
0°. The isotopic ratio decreases with increasing 6, be- 
cause for the same solid angle as viewed from the star, 
the 26 Al atoms ejected at backward angles are mixed 
with a higher mass of shocked CS gas. The mass con- 
tained in the shell's forward hemisphere (i.e., at 8 < 90°) 
is nevertheless significant: it amounts to 23 Mq for the 
25 Mq star and is between 200 and 250 Mq for the three 
other stars. 

Hydrodynamic simulations of stellar wind bow shocks 
have shown that the steady-state solution of IWilkinl 
(1996) provides a good description of the time-averaged 
shape of the bow shock shell; although a bow 
shock is neither smooth nor steady (|Raga et al.l 119971 ; 
iBlondin fc Koerwer1ll998D . The shell is subject to peri- 
odic oscillations in and out with respect to the equilib- 
rium position, which has been interpreted as resulting 
from the nonlinear thin shell instability (NTSI; I Vishniaa 
Il994f h This instability is also known to be an efficient dy- 
namical focusing mechanism for large-scale gas stream, 
resulting in the buildup of dense cores. Thus, hydro- 
dynamic simulations of the NTSI have shown that the 
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t c (yr) 

Fig. 3. — 26 A1/ 27 A1 ratio in a prestellar core formed in the 
shocked-gas layer of a SN remnant resulting from the explosion 
of a runaway WR star, as a function of the time of core collapse 
after the explosion. Same as for Figure 2, V* = 20 km s _1 and 
rtjj = 100 cm' 3 . The hatched areas reflect a factor of 2 of uncer- 
tainty in the calculated isotopic ratio. The dotted horizontal line 
shows the canonical ratio measured in CAIs. The results for the 
25 and 60 Mq stars are almost identical. 

density contrast in the shell can reach 10 2 to 10 4 , depend - 
ing on the gas cooling efficiency (e.g.|Hueckstaed3[2003). 
The high-density seeds thus generated are l ikely sites of 
further star formation (He itsch et al.|[2008l ). However, 
the required gravitational collapse of these dense cores 
is probably not possible as long as the shell is exposed 
to the intense photoio nizing radiation from the nearby 
massive star (see, e.g. JRaga et al.iri997t ). 

At the end of the WR stellar phase, the SN outburst 
will expel the bow shock material to large distances. 
About 10 4 years after the explosion, radiative cooling of 
the shock-heated gas will become important and the SN 
remnant will enter the pressure-driven snowplow phase. 
The transition from the adiabatic to the radiative phase 
in SN remnants is accompanied by the development of 
dynamical instabilities, that can further increase th e 
mass of pre-existing gas clumps (|Blondin et al.l [l998h . 
The timescale for collapse of a dense core embedded in a 
shocke d gas layer is govern ed by the gr avitational insta- 
bility (Heitsch et al.l l2008h and reads (jWhitworth et al.1 
Il99"l 



2c„ 



(2xl0 6 yr) 



10 22 cm -2 



(7) 



where S coro is the H column density of the core, c s ks 
0.5 km s -1 is the local sound speed, and G is the gravi- 
tational constant. 



The 26 A1/ 27 A1 ratio in such a prestellar core can be 
estimated to be 



/ 26 A1\ 



26 A1 \ 



■ exp(-i c /r 26 ) . (8) 



This ratio is shown in Figure 3 as a function of t c . We 
adopted for ( 26 A1/ 27 A1) S and E s the values at 9 = 45° 
(see Figure 2). The error in ( 26 Al/ 27 Al) coro shown in 
Figure 3 is intended to account for various uncertainties 
in the model parameters, e.g., in 9, A r 2 6, and c s . We see 
that ( 26 Al/ 27 Al) core increases almost linearly up to t c 



10 yr, thus reflecting that E c 



oc t~ 



(Equation 0) 



The predicted 26 A1/ 27 A1 ratios are consistent with the 
canonical value measured in CAIs for a large interval of 
t c (as the delay for the CAI formation after collapse of 
the presolar nebula is <C r 26 , it can safely be neglected). 
The exponential decay of ( 26 Al/ 27 Al) coro for t c > 10 6 yr 
is due to the 26 Al radioactivity. 

Inserting Equations (J5J) and © into Equation © one 
can see that ( 26 Al/ 27 Al) core scales as n H X V^. Thus, 
the solar system formation could have been triggered by 
the explosion of a runaway WR star propagating into 
atomic interstellar gas (jih ~ 1 cm -3 ), but only if the 
star's velocity was > 100 km s _1 . 

In the proposed scenario, the incorporation of 26 Al- 
rich stellar ejecta into interstellar gas is due to various 
dynamical instabilities operating in both the bow shock 
and the SN remnant shells. The associated turbulence 
is expected to homogenize the mixing at all scales, re- 
gardless of the carrier phase of 26 Al (gas or dust). This 
is con sisten t with the Mg isotopic data of iThrane et al.l 
(|2006fl and iVilleneuve et all (|2009H which suggest that 
26 Al was homogeneously distributed in the early solar 
s ystem. 

lArnould et all (|2006l ) showed that WR star nucleosyn- 
thesis can produce 26 Al, 36 C1, and 41 Ca at levels compat- 
ible with the meteoritic measurements, provided that the 
delay before the incorporation of these SLRs into CAIs 
was ~l-3 x 10 5 yr. But using the yields for a 60 Mq star 
given by these authors, we obtain a 36 CI abundance well 
below th e value reported in CAIs, as also found previ- 
ously by iGaidos et all (|2009f ). On the other hand, the 
present work shows that both 26 Al and 41 Ca abundances 
in meteorites can result from the contamination of the 
presolar molecular core by material ejected from a run- 
away WR star, whose explosion as a SN triggered the 
formation of the solar system. 

We thank Jean-Pierre Thibaud, Fai'rouz Hammache, 
and Pierre Roussel for their critical reading of the 
manuscript. 
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